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Abstract 

We performed time lag analysis on the X-ray light curves of Centaurus A (Cen A) obtained 
by the Gas Slit Camera (GSC) aboard the Monitor of All-sky X-ray Image {MAXI) in three 
energy bands (2-4 keV, 4-10 keV, and 10-20 keV). We discovered a soft X-ray lag relative to 
higher energies (soft lag) on a timescale of days in a flaring episode by employing the discrete 
correlation function (DCF) and the z-transformed discrete correlation function (ZDCF) method. 
In the episode, a peak and a centroid in the DCF and the ZDCF was observed at a soft lag of 
5 days in 2-4 keV versus 4-10 keV and in 4-10 keV versus 10-20 keV, and ~ 10 days in 2-4 
keV versus 10-20 keV. We found it difficult to explain the observed X-ray variation by a single 
energy injection with the one-zone synchrotron self-Compton (SSC) model, in which the soft 
lags in these three energy bands reflect the different cooling times of the relativistic electrons, 
by assuming the magnetic field and minimum Lorentz factor estimated from a broad-band 
SED. Alternatively, if the phenomenon is interpreted as cooling of Comptonizing electrons in 
a corona covering the accretion disk, the temperature of the corona producing the variable 
X-rays should be ^10 keV for reconciliation with the soft lag in the energy range of 2-20 keV. 


1 



Keywords: Galaxies: active — Galaxies: Seyfert — X-rays: individuals (CentaurusA) 


1 Introduction 

Centaurus A (Cen A, NGC 5128), at a distance of 3.8 Mpc (Harris et al. 2010), is the nearest aetive 
galaetie nuelei (AGN) from the Milky Way galaxy and is known as the brightest Seyfert galaxy in 
the X-ray sky (2.87 x 10“^° erg cm“^ in the 4-10 keV band; Hiroi et al. 2011, and 7.48 x 
IQ-io in the 14-195 keV band; Tueller et al. 2008). In addition, this objeet has often 

been eonsidered as a jet-dominated AGN seen from the side ealled “misdireeted BL Lae object” 
(Morganti et al. 1992; Chiaberge et al. 2001) with its jet inelined to the viewing angle at ^ 50° 
(Tingay et al. 1998). Owing to its brightness and proximity, this objeet has been an ideal target to 
study the geometry and the emission mechanism around a supermassive black hole. However, the 
origin of the X-ray radiation from Cen A has not been easily identified from its energy speetrum. The 
observed X-ray speetrum has no significantly different features from typical Seyfert galaxies, whieh 
radiate X-rays mainly from the inner disk region with Comptonization in hot plasmas loeated above or 
around the aeeretion disk {e.g. Rothsehild et al. 1999; Beckmann et al. 2011; Fukazawa et al. 2011). 
On the other hand, the high energy tail of the X-ray speetrum appears to be smoothly eonneeted to 
the MeV emission whieh is believed to be generated from a jet, and the speetral energy distribution 
(SED) in the energy range between radio and GeV 7 -ray ean be modeled sueeessfully by the one-zone 
SSC model or other varieties of jet emission models {e.g. Abdo et al. 2010; Roustazadeh & Bottcher 
2011; Petropoulou et al. 2013), though a eombination of a thermal inverse Compton eomponent and 
a non-thermal jet emission also describe the SED equally well (Abdo et al. 2010). 

When the speetral analysis does not provide a unique emission meehanism or physical pa¬ 
rameters of AGNs, time variability analysis, espeeially time lag analysis, may be useful to obtain 
additional information to resolve the issue. Kataoka et al. (2000) deteeted soft lags {i.e. the variation 
in a lower energy band lags respeet to that of a higher energy band) on a timeseale of lO^-lO^ sec 
in an X-ray flare from the TeV blazar PKS 2155-304 and estimated the magnetic field strength at the 
X-ray emitting region in the jet with the hypothesis that the delays are caused by different synehrotron 
lifetimes of the relativistie eleetrons. Similar diseussion was presented by Tanihata et al. (2001) for 
other blazar objeets: Mrk 421 and Mrk 501. Sriram et al. (2009), eonversely, found hard lags (vice 
versa for soft lag) on a timeseale of ~ 10^-10^ sec in the Seyfert 1 galaxy NGC 4593 using the XMM- 
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Newtonl^VlC. The authors interpreted the hard lags as a timeseale of thermal Compton seattering in 
a eompaet eentral eorona and derived its size as about five Sehwarzsehild radii {Rs) for the blaek hole 
mass 7 x Recently, De Marco et al. (2013), based on the study of 32 radio quiet AGNs, dis¬ 

covered the relation between black hole masses (Mbh) and soft X-ray lags (r ~ 10-600 sec) in high- 
frequency variation (z/ 0.07-4 x 10“^ Hz); log|r| = 1.98[±0.08] -f 0.59[±0.11]log(MBH/10'^MQ). 

For Cen A, time lag analysis in the X-ray band has scarcely been applied as yet, presumably 
because of remarkably stable spectral (or flux) behavior as reported by Rothschild et al. (2011) based 
on 12.5-year observation by the Rossi X-ray Timing Explorer combined with pointed observation 
with inadequate time coverage or continuous monitoring with insufficient sensitivity in the X-ray 
band. Monitor of All-sky X-ray Image (MAXI; Matsuoka et al. 2009), launched in August 2009, has 
a high sky coverage of 95% and a X-ray sensitivity to detect above 15 mCrab in the 4-10 keV band 
with a cadence of ~90 min (Sugizaki et al. 2011), and thus allowing us to investigate a long-term 
activity for AGNs with a sufficient time resolution. 

In this paper, we study X-ray variation properties in Centaurus A. The lags in long-look flux 
variations are investigated using light curve data sets obtained by the MAXI/GSC in MJD = 55058- 
56800 in §4.1. Then we focus on a flare-1 ike enhancement event observed during MJD = 55075- 
55160 in §4.2. For the lag analyses, we employed the continuously evaluated discrete correlation 
function (CEDCF) based on the discrete correlation function (DCF; Edelson & Krolik, 1988), and 
the z-transformed discrete correlation function (ZDCE; Alexander, 1997). Einally, we discuss the 
interpretations of this phenomenon and the origin of X-ray emission of Cen A in §5. 

2 Observations and Light Curves 

In August 2009, MAXI was launched and attached to the Japanese experiment module “Kibo” on 
the International Space Station (ISS). Since then, the gas slit camera (GSC; Mihara et al. 2011) and 
solid-state slit camera (SSC; Tomida et al. 2011) aboard MAXI has been producing an all sky X-ray 
image at every orbital period of ISS (92 minutes). The SSC employs X-ray CCD arrays and covers 
the energy range from 0.5 keV to 12 keV, while, the GSC employs gas proportional counters and 
covers the energy range from 2 keV to 20 keV. Eigure 1 shows the MAXIIGSC light curve during 
MJD=55058-56800 in the three energy bands obtained by the MAXI on-demand analysis (ver 2.0) 
implemented at the MAXI public web site‘, where calibration described in Nakahira et al. (2012) is 
applied. These light curve are binned by five days. The flux histories of Cen A reveal variations of a 
factor of ~3-5 on a timescale of the order of several weeks or months. In this analysis, we used only 


‘ http://maxi.riken.jp/top/ 
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the bin meeting the eondition that the data whose error size is smaller than e* + Sag., where the 
parameter dj is the mean of e*, and is the standard deviation of Cj. The purpose of this treatment is 
removing the low quality data resulting from a relatively insuffieient number of seans in the five day 
interval. 5-7% of the data points were removed by the treatment deseribed above. It should be noted 
that the sensitivity of the MAXI/GSC has deteriorated over time due to damages in the deteetors and 
ehanges of the operation modes. For more information, see Mihara et al. (2014). 

3 Method of Time Lag Analysis 

To evaluate time lags in X-ray variation between two different energy bands, we employed two teeh- 
niques based on the diserete eorrelation funetion, ealled the eontinuously evaluated diserete eorrela- 
tion funetion and the z-transformed diserete eorrelation funetion. 


3.1 Discrete correlation function 


The DCF was introdueed by Edelson and Krolik (1988). To analyze the irregularly sampled astro- 
nomieal data, DCF has been used widely to uneover eorrelations or lags between flux variations in 
two different energy bands, espeeially for AGNs. This funetion provides the eorrelation eoeffieient 
between two unevenly sampled time-series data as a funetion of lag (riag), and the value riag cor¬ 
responding to the peak or the centroid of DCF{tu^), is presumed to be the value of the lag in the 
variations in two time series. First, we calculate the set of un-binned discrete correlation functions 
(UDCFs) for all measured pairs (ai,bj) defined as follows: 


UDCFij = 


a bj —b 
— X -, 


( 1 ) 

where d and b are the means, and a a and ab are the standard deviation of data sets a and b. Each 
UDCFij is associated with the pairwise lag AUj = tj — U. Next, we average the M pairs for which 


over Tlag ^r'lag/2 ^ Atij Tlag T Ariag/2, 


DCF{T^) = -J^UDCF,j. 




One a error of each bin is defined by the following expression: 

1 1/2 


f^(Aag) = 


M-1 


j:{UDCF,,-DCF{n,,)y 




( 2 ) 

( 3 ) 


3.2 Continuously evaluated discrete correlation function 

The continuously evaluated discrete correlation function (CEDCF) was originally introduced by 
Goicoechea et al. (1998). The CEDCE is an oversampled DCE, meaning that the bin-to-bin inter- 







val is smaller than the bin width, and thus this method allows us to determine more precise time lags 


from the position of the peak or centroid than via the standard DCF method. Details of advantages of 


this method are described in Gil-Merino et al. (2002). Here, we used the centroid instead of the peak 
lag, because the centroid is less sensitive to a particular binning used and statistical fluctuation. In 
order to estimate the significance and uncertainties of the time lags, we built a cross-correlation cen¬ 
troid distribution (CCCD) by running Monte-Carlo simulations known as “flux redistribution/random 
subset selection” (FR/RSS) described in detail by Peterson et al. (1998). For every simulation, in 
order to account for the effects of uncertainties in the measured flux, we altered each measurement 
value to a random Gaussian deviate by assuming that flux uncertainties are Gaussian distributed about 
measured value with uncertainty (FR), and in order to account for the effects of sampling in time, we 
select random N points from N real data points regardless of whether they have been previously se¬ 
lected or not (RSS; known as “bootstrap method”, and this procedure reduces the number of points 
in each light curve by a factor of ~ 1/e). This method yields fairly conservative uncertainties in the 
obtained lags, i.e., the real uncertainty may be somewhat smaller (Peterson et al. 1998). Here, we 
calculated the centroid of the CEDCF (r^^™*^^) as {TjiTiCEDCFi)/{T^iCEDCFi) using the points 
with correlation coefficients larger than the half of the peak {CEDCFi > 0.5 x CEDCFpeak)- In 
each simulation, the centroid of a CEDCF is determined and recorded to build the CCCD, and then 
we can obtain Tmedian and iAres directly from the constructed CCCD respectively, where ± Ares are 
defined as that 15.87% of the realizations give values below Tmedian — Ares and above Tmedian + Ares 
{i.e., 68.27% of yielded centroids are contained between Tmedian — Ares and Tmedian + Ares). Thus, 
the values iAres correspond to la errors for a normal distribution. 

3.3 Z-transformed discrete correlation function 

Alexander (1997) proposed the z-transformed discrete correlation function (ZDCF) to improve the 
CCF analysis of sparse and unevenly sampled light curves. Several biases of the conventional DCF 
method can be corrected by the ZDCF technique by using Fisher’s z-transform and equal population 
binning rather than equal time interval Ariag. If we have n pairs in a given time lag riag, the correlation 
coefficient at this lag r(riag) is estimated by 



( 4 ) 


where a,h and aa,af, are the means and the standard deviations of data sets a and b, respectively. 
The sampling distribution of r is highly skewed {i.e, the distribution is far from normal distribution), 
and therefore calculating its sampling error by simple standard variance of r might be very inaccu¬ 
rate. Assuming a and b are drawn from the bivariate normal distribution, we can transform r into an 
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approximately normally distributed random variable by Fisher’s z: 



( 5 ) 


The mean and the varianee of 2 ; are approximately equal to 



n-l[ 2(n-l) 6(n-l)2 


5 + p2 ^ Il + 2p2 + 3p4 ^ 


( 6 ) 


(7) 


These relations yield the interval corresponding to the ± la error interval which is represented as 


5r± = I tanh(z ± s^) — p\■ 


( 8 ) 


The binning is performed with a minimum lag width of e and a minimum UDCF number of rimin 
for each bin, where is at least 11 points for a meaningful statistical interpretation. Also, the 
interdependent pairs are discarded in each bin; namely a new pair whose a or 6 points have previously 
appeared in that bin is not gathered anymore. The lag time associated with a bin is defined as the mean 
value in the bin, and its upper and lower errors (Ar±) are defined as the intervals that contain 34.1% 
of the points in the bin above and below f, respectively. It is more robust method than conventional 
DCF method when the light curves are very sparsely and irregularly sampled. The ZDCF time lag is 
characterized by the centroid which can be calculated through the same procedure as that of 

the CEDCF. 

4 Analysis and Result 

4.1 Analysis of Complete Time Series 

First, we considered the 4.5 years light curve data sets. To search for possible time lags and correlation 
between each fluctuation in the three energy bands, we calculated the CEDCFs(riag) and ZDCEs(riag) 
for all energy combinations (2-4 keV vs. 4-10 keV, 4-10 keV vs. 10-20 keV, and 2-4 keV vs. 10-20 
keV) for the range of lags between —200 days and +200 days, where the CEDCEs are calculated with 
the following values: 5 days as the bin-to-bin interval and 15 days as the bin width Ar^ag (see §3.2). 

Results are displayed in Eigure 2. We present the CEDCEs and the error ranges associated with 
them as the gray shadowed areas, and the ZDCEs as black circles with error bars. Here the CEDCEs 
and the ZDCEs show almost the same behavior. It is to be noted that all of the maximum values of 
the correlation coefficients calculated with the both methods are located at Lag > 0. Again, in this 
paper, a positive time lag indicates that the variation of soft X-rays is delayed with respect to that 
of hard X-rays. This result, therefore, suggests that a “soft lag” on timescales of several days exists 
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in flux variations. We will examine the signifleanee of the soft lags quantitatively below. A similar 
lag on timeseales of a day in the X-ray band has been reported in NGC 4151 by Caballero-Gareia 
et al. (2012). They found the asymmetrie shape of DCF between light eurves in the 20-50 keV and 
50-100 keV obtained by Swi///BAT, i.e. DCF{t\^^ > 10 days) is larger than DCF{Tia,^ < —10 days), 
suggesting that the hard band lead the soft band variations on timeseales longer than ~ 10 days. 

To evaluate the statistieal uneertainties of the soft lags, we performed Monte-Carlo simulations 
utilizing the FR/RSS method. For eonstrueting CCCDs, we ealeulated CEDCFs for 1000 light eurves 
produeed by Monte-Carlo simulations with a bin size of 10 days. In order to minimize the systematie 
effeet of bin size, we performed similar ealeulations with 15, 20, and 25 day bins, and eonsequently 
a total of 4000 simulations are performed to obtain a distribution of eentroids. In the same manner 
as those for the CEDCF, we also performed simulations for smoothed ZDCF ealeulated by using 
running average of 10, 15, 20, and 25 days, where we seleeted rimin = 11 and e = 5 days for eaeh 
simulated ZDCF in this analysis. The distributions of the eentroid of CEDCFs and smoothed ZDCFs 
are presented in Figure 3. Solid lines show the CCCDs derived through the CEDCF method, and 
dotted lines show the CCCDs built with the ZDCF method. On the basis of these distributions, we 
ealeulated the 1 a error ranges of the soft lags. The results are eonsistent between two methods. Table 
1 shows the magnitudes of the soft lags and assoeiated la error estimated. The first eolumn gives the 
energy bands used for the ealeulation, and the seeond and fourth give the possible lag obtained by the 
CCCD of a CEDCF and a ZDCF respeetively. Based on the CCCD, we also show the 

probability P<o that the true lag is less than zero in the table. Here, the realizations of P<o gives the 
probability of the simulated eentroids lie in the negative day region, where soft lag is represented as a 
positive value. In spite of the eentroids with highest possibility in 4-10 keV vs. 10-20 keV and 2-4 
keV vs. 10-20 keV are evaluated to be positive values in both methods, their error ranges inelude the 
zero, i.e. no lag, and the eonfldenee levels of the soft lags are less than or eomparable to 82.8%. The 
soft lags averaged over the eomplete time series, therefore, are not highly eonfldent. 

Next, in order to elarify the possible time lags and eorrelation between eaeh fluetuation in 
the three energy bands as a funetion of time, we ealeulated the running CEDCF(tk,Tiag) between 
separated data sets (a*, bj) eontained in the time width tk — At/2 < t < tk + At/2 (k = 0,1,2 • ■ •) 
by employing a reetangular window funetion. The resulting running CEDCF(ffc, Tiag) is depieted in 
Figure 4. Here, we set the parameter At as 200 days, Ariag as 20 days, and the distanee between 
eenters of bins (tk+i — tk) as 5 day. The horizontal axis eorresponds to the eenter of the window 
(tk), and the gray seale denotes the eorrelation eoeffleient. To assist eomparison, the light eurve of 
the 2-20 keV is also displayed on the upper side. Until MID ~ 55800 the CEDCFs have relatively 
high eorrelation eoeffleient (represented by dark gray) almost eontinuously in all CEDCFs. Here, the 

7 


Energy band (keV) (days) P^o^days) (day) P^f/^days) 


CEDCF 


^ZDCF , 


^ZDCF 


2-4 vs 4-10 

0 1+'^® 
*^•^-4.5 

0.484 

+ + 

+ 1 

CO 

o 

1 

0.536 

4-10 vs 10-20 

1 q+4.8 
-*-•'^-4.5 

0.314 

1 1+3.6 
-*-•-*--3.0 

0.359 

2-A vs 10-20 

Y 9+7.6 
''^-7.6 

0.172 

N 

CO 00 

+ 1 

O 

0.253 


Table 1. The estimated soft lags using FR/RSS for CEDCF and ZDCF methods. In this analysis, CEDCFs are calculated with one day 
intervals and nmin is confined to 11 for calculating ZDCFs. The details of these methods are written in the text. 


peak of the correlation coefficient stay roughly constant in all cases. The lag is ~ 0 days for 2-4 
keV vs. 4-10 keV, while ~ 10 days for the others. These lags are consistent with the results shown in 
Figure 2. Especially, the portion between 55075-55150, which contains the highest flux enhancement 
event, yielded rahter clear peaks of the CEDCF in all energy combinations. On the other hand, there 
are not prominent features in CEDCFs during the period MJD ~ 55800 to MJD ~ 56250, probably 
due to insignificant flux fluctuations in each time window of 200 day. For MJD > 56250, while the 
coefficients are lager than the previous period as a whole, peaks of DCFs are broader and centroids are 
rather ambiguous. Although the changes of the shape and the feature in the CEDCFs might reflect the 
transition of X-ray emission mechanism in Cen A, the episodes in the later epochs are more difficult 
to study due to reduced sensitivity of the MAXI/GSC. In this parer, we focus on the first flare-like flux 
enhancement episode, where the correlations are particularly significant. 

4.2 Analysis of the flux enhancement event 
4.2.1 Light Curves 

On the lag analysis described above, we found that the strong surge of X-ray flux observed during 
55075-55150 contributed to the soft lag suggested in §4.1. For studying this event more reliably, we 
remade the light curve by using the maximum-likelihood image fitting method as described in Hiroi 
et al. (2013). This method allows us to suppress the back ground fluctuation, and thus make a high 
quality light curve. We employed the following procedure with the data selection criteria and the 
background model described in Hiroi et al. (2013). We excluded the event data taken at high latitude 
regions (|latitude] > 40°), where the charged particle flux is high, and those detected in near both 
ends of each proportional counter at a photon incident angle \(j)\ > 38° (for the definition of 0, see 
Mihara et al. 2011) where the background rejection efficiency is lower. While this refined background 
model is available only in a limited to a portion of time in the MAXI observation in energy range of the 
2-4 keV and the 10-20 keV at present, this method successfully reduces the flux errors attributable to 
background fluctuations. The remade light curves are not significantly different from the one obtained 
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through the MAXI on-demand proeess in this epoch except for uncertainties in each bin. 

Figure 5 shows the light curve of Cen A with a bin size of one day during MJD=55075-55160 
in the three energy bands, beneath which we show the hardness ratio (the 4-10 keV divided by the 
2-4 keV). Calculated fluxes binned in 1 day are represented in black filled circles. We can see a 
flux enhancement peaking around MJD ~ 55090 at a flux ~ 60 mCrab in the 4-10 keV band. It is 
the brightest flux recoded by the MAXIIGSC so far for this object. The relative amplitudes of flux 
enhancements (a factor of ~ 2) are similar in the three energy bands. On the other hand, the temporal 
variability looks somewhat different among the three. The flux decay time in the 2-4 keV seems 
longer than those in the other energy bands, while the timescale of variation in the 10-20 keV seems 
shorter than those in the other energy bands. As for the hardness ratio in figure 5, we can see the 
decreasing trend toward the peak of the flux enhancement event from MJD = 55075 to MJD = 55100, 
and then it stays almost at a constant. 


4.2.2 Time lag Analysis 

Using the method of the time lag analysis based on the DCF described in §3, we investigated the cor¬ 
relations and time lags of the X-ray fluctuations among three energy bands during the flare observed 
at MJD ~ 55090. The CEDCFs and the ZDCFs between the 2-4 keV and 4-10 keV, the 4-10 keV 
and 10-20 keV, and the 2-4 keV and 10-20 keV are shown in Figure 6. The CEDCFs calculated with 
a bin size of 2 days and a distance between centers of bins of 1 day are shown by shadowed areas, 
and the data points and the error bars plotted on the shadow areas represent the ZDCFs. One can see 
in Figure 6 that the calculated cross-correlation functions through the two methods give almost the 
same results. The centroids of CEDCFs are located at -1-7.3 days, -1-6.7 days, and -1-11.4 days, while 
those of ZDCFs are found at -1-3.1 days, -1-1.4 days, and -1-7.8 days for the three pairs, respectively. 
These results indicate that the existence of soft lags of ~5-10 days in X-ray flux variation from Cen 
A. Additionally, the shift of the DCF centroid to a positive direction {i.e., soft lag) are caused by the 
overall shape of the DCF, in other words, it is attributable to the trend on a timescale of ~ 20-30 days 
(estimated from the HWHM of the DCFs), rather than to the small fluctuations on a timescale of ~ 
3-5 days. The soft lags, therefore, arise not from the rapid fluctuations especially found in the light 
curve of the 10-20 keV but rather from the more slow variations on a similar timescale as the HWHM 
of DCFs. Correspondingly, even if we use the light curve with the bin size of three day for evaluating 
DCFs, the features described above are reproducible. 

The statistical uncertainties of the soft lags are derived through the same procedure as that used 
in§4.1, i.e. the FR/RSS method, by changing the bin width from 2 day to 5 day with a step of 1 day. 
For ZDCF, we set the two parameters as nmin = 11 and e = 1 day. The evaluated results from 4000 
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Energy band (keV) 

(days) 

pCEDCF 
■^<0 (days) 

(day) 

pZDCF 
-^<0 (days) 

2-4 vs 4-10 

4-10 vs 10-20 

2-A vs 10-20 

6.7tll 

6.3111 

12-2119 

0.072 

0.104 

0.008 

4-9lli 

3-4llg 

9-011*9 

0.089 

0.202 

0.011 


Table 2. The estimated soft lags during MJD = 55075-55160 using FR/RSS for CEDCF and ZDCF methods. In this analysis, CEDCFs 
are calculated with one day intervals and is confined to 11 for calculating ZDCFs. The details of these methods are written in the 
text. 


Monte-Calro simulations are summarized in Table 2. Confidenee levels of soft lags are signifieantly 
improved eompared with the ease of the whole light eurves, and reaeh 99% in the 2-4 keV vs. 10- 
20 keV in this period. Again, the FR/RSS method yields a somewhat broadened distribution of the 
eentroid of DCF compared to the true situation. Thus, these results, especially the case of the DCF 
between the 2-4 keV and 10-20 keV, suggest that the day-scale soft lags exists in X-ray emission 
from Cen A at the confidence level of potentially > 99% in flux variations on a timescale longer than 
20 days. 


5 Discussion 

We have found a soft lag of ~ 5-10 days in Cen A in the X-ray flux variations on a timescales longer 
than ~ 20 days between the 2-4 keV and 10-20 keV bands with a confidence level of 99% in the flux 
enhancement event at MJD ~ 55090. It could be argued that the high significance of the soft lag in the 
first flare episode should be discounted because the significance in the later episode are not so high. 
We note, however, that the first episode had the highest peak flux in the 4-10 keV (and also 2-20 keV), 
and that the sensitivity of the MAXI/GSC was best at that epoch. While it is difficult to demonstrate 
that a soft lag is always present in every flux enhancement events, the boxcar CEDCF for the later 
episode in Figure 4 shows a similar feature to that of the first event, and suggestive of the existence 
of a soft lag. In the following sections, we discuss if the observed features described above are 
accountable by the X-ray emission or reprocess scenarios for AGNs; namely, the synchrotron cooling 
of electrons in the one-zone synchrotron self-Compton process in the jet, cooling of electrons by the 
Compton up-scattering of the disk photons in the corona, and X-ray reverberation on the accretion 
disk or a warm absorber. 
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5.1 One-zone synchrotron self-Compton 

Soft lags in AGNs are sometimes explained by the radiative cooling of relativistic electrons through 
the one-zone synchrotron self-Compton (SSC) process (e.g. Kataoka et al. 2000). In this model, the 
power produced by a single electron in the synchrotron radiation and the inverse Compton scattering 
of a single electron is given by (Rybicki & Lightman 1979) 


^ = ^[(7 - l)meC^] = ^CCTt( 7 ^ - 1){Ub + Cph), 


(9) 


where rrie is the mass of an electron, ctt is the Thomson cross section, and U-q and f/ph are the energy 
densities of magnetic field (B'^/Sir) and the seed photon field, respectively. The cooling timescale, 
which sets the characteristic timescale for the kinetic energy loss of electrons by the synchrotron 
radiation and the inverse Compton scattering, can be expressed as follows: 


fcool 


(7- l)meC^ 


SrrteC 


( 10 ) 


4ccrT(7^ - l)(f/'B + C'ph)/3 4 (Tt( 7 + !)(% + C'ph) ’ 

In the jet comoving frame, the time lag (tiag) of emission at 72 with respect to emission at 71 is the 
time it takes for electrons to lose energy equal to (71 — 72 )meC^. Thus, flag in the jet comoving frame 
is given by 

3nipC / 1 1 \ 


( 11 ) 


( 12 ) 


® 4(Jt(Cb + Cph) V 72 + 1 7i + 1 y 

In this paper, we consider the upper limit of the cooling time of an electron: 

SrrteC 1 
® 4(Tt(7b 7min + 1 ’ 

where 7min is the minimum Lorentz factor of electrons in the emission region. Abdo et al. (2010) 
attempted to explain the multi-wavelength SED of Cen A with the one zone SSC model and yielded 
the magnetic field of the emitting region B = 6.2 G, a minimum electron Lorentz factor 7 jnin = 300, 
and the doppler factor <5 = 1. Assigning these values to Eq.(12) and considering that the cooling time 
in the observer’s frame is not different in the comoving frame under the situation of <5 = 1 , we obtain 
< 1.1 day as a conservative upper limit of radiative cooling time of the electrons in the observer’s 
frame. The upper limit seems only marginally consistent with the observed day-scale soft lags within 
flux variations on a timescale ^ 20 days, if the observed flux enhancement is caused by a single energy 
injection. One way to produce a variation timescale longer than an intrinsic electron cooling time is 
the emission region with a light-crossing time longer than the fcooi- In this case, however, the shape 
of light curves expected be close to symmetric as suggested in Kataoka et al. (2000), in contradiction 
to the X-ray light curves shown in figure 5. Additionally, an emission region size of ~20 light days 
corresponding to ~850 Schwarzschild radius (Rs) for the black hole mass of Mbh = 2.0 x 10 ®Mq 
(M arconi et al. 2001) is required to construct such a condition. It seems to be too large for the emitting 
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region size ~ 3 x 10^® em (eorresponding to ~ 50i?s) estimated by fitting the broad-band SED (Abdo 
et al. 2010). If the enhaneement of flux intensity with soft lags was triggered by a single energy 
injeetion, therefore, a weaker B field or/and a small 7 min are necessary at the period. 

Alternatively, multiple injections may be able to produce asymmetry in a CEDCF and a soft 
lag longer than the fcooi at the same time if the time resolution and photon statistics in the present 
data is not sufficient to resolve individual flares. Eor example, hard-to-soft evolution in hardness ratio 
is often observed in gamma-ray bursts (GRBs) with multiple pulses, both within a single pulse and 
over series of pulses in a burst {e.g. Eishman & Meegan 1995, Eenimore et al. 1995). If multiple 
energy injections in a homogeneous region {e.g. series of internal shocks caused by repeated ejection 
of matter), or coincidentally discrete single energy injections in several regions, occurs in Cen A 
and undergoes hard-to-soft evolution, a soft lag could be observed. In fact, a decreasing trend in the 
hardness ratio (i.e the softening trend) is suggested around the peak of the flux in Figure 5. SSC 
process with a multiple injection scenario is, therefore, not rejected by our results. 


5.2 Up-scattered black body radiation 

The X-ray emission in a Seyfert galaxy is commonly thought to be produced by inverse Compton 
scattering in a hot corona surrounding the inner parts of the accretion disk. Here we discuss the 
reasonable situation that can yield the day-scale soft lag in variations on a timescale ~ 20 days on the 
basis of this assumption. 

In the Comptonization process, the variation in flux may be caused either by the seed photon 
or by the Comptonizing plasma. The former case, naturally leads to hard lag associated with the 
propagation of mass accretion fluctuations. A soft lag, on the contrary, may be caused by the cooling 
of Comptonizing plasma. However, the hot corona implied by the main X-ray power-law continuum 
which extends to hundreds of keV, seen in INTEGRAL data of Cen A (Beckmann et al. 2011, Burke 
et al. 2014), is not relevant to produce variations with soft lag in the soft X-ray region (2-20 keV). To 
produce the observed soft lag, the plasma temperature (or a spectral break energy) needs to be around 
~ 10 keV. However, we think it difficult to expect the temperature of the hot corona responsible 
for the X-ray continuum to decrease from hundreds of keV to ^ 10 keV considering the stability 
in the X-ray band as observed by RXTE (Rothschild et al. 2011). Such a low-temperature corona, 
alternatively, is possibly present in addition to the hot corona mentioned above, and its cooling may 
cause the observed soft lag. This idea has been applied to explain the Comptonized component with a 
low break energy found in a Seyfert galaxy Mrk 509 (Noda et al. 2011). They found that a soft excess 
component detectable at ^ 3 keV varies independently of the main X-ray continuum, and argued that 


12 


this soft component with kTe = 0.49 keV is a evidence for “multi-zone Comptonization” similar to the 
case of the black-hole binary Cygnus X-1. In Mrk 509, the soft excess component was stable during 
a few days but varied on timescales of a few weeks. If Cen A has similar multi-zone Comptonization 
component and the soft component has the plasma temperature of several keV, observed soft lag 
may be explained by cooling of the Comptonizing plasma. The nuclear emission of Cen A is not 
detectable in optical or UV spectrum, due to a high column density (A^h ~ 1-6 x 10^^ atms/cm^; 
Rivers et al. 2011). Therefore, we do not have good knowledge of the seed photon from the inner 
accretion disk, which influence the cooling time. The cooling time of the Comptoniazation plasma in 
the non-relativistic limit, however, is suggested to correspond to the power spectrum density (PSD) 
break in the X-ray band of Seyfert galaxies (Ishibashi and Courvoisier, 2012). For Cen A, Rothschild 
et al. (2011) measured the PSD break in 1.5-12 keV band at Tb = 18.3^g®7^ days. It is consistent 
with the timescale of the variation in our results. These argument supports our interpretation that the 
observed soft lag is due to a soft component of multi-zone Comptonization. 


5.3 Reverberation on the Accretion Disk, Warm Absorbers, the BLR, or the Torus 

Previous observations have revealed that X-ray spectra of AGNs often show a soft X-ray excess below 
~ 1 keV, superposing above a main power law component. This phenomenon is usually interpreted 
as the reflection paradigm where a corona (power-law component) irradiates the warm absorber lying 
around 100i?s (corresponding to ~ 2.3 light days for Cen A) from the black hole {e.g. Miller et al. 
2008 ), or the inner region of the accretion disk within 20i?s (corresponding to ^ 0.46 light days for 
Cen A) where relativistic blurring is significant {e.g. Cackett et al. 2013). In the both hypotheses, the 
soft excess is naturally explained as the contribution of Fe L lines in these origin. 

Recently, De Marco et al. (2013) reported lags of soft excess with respect to the X-ray con¬ 
tinuum in over a dozen AGNs, and explained that the delayed soft excess emission originates in the 
innermost regions of the accretion disc. Our data, however, starts at 2 keV and do not contain the en¬ 
ergy region containing the Fe L lines. This model, therefore, predict a hard lag in contradiction to our 
results on the contrary. With reflection, the variation of the 2-4 keV continuum flux that represents 
the direct component lead the variation in the 4-10 keV range containing the Fe K line and the 10-20 
keV band containing the Compton hump. In fact, Zoghbi et al. (2014) reported a time delay between 
energies > 10 keV and the continuum in MCG-5-23-16. They claimed that the these lags are most 
likely due to reverberation in the reflection Compton hump in a manner similar to the response of the 
relativistic Fe K line. Besides, a reflection component from the broad line region (BLR) or the torus 
located at a range of a mili-parsec to a parsec (a few light days to a few light years; e.g. Peterson 
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1993, Krolik and Begelman 1988) from the center would be detected in the X-ray spectrum as a hard 
lagging component. Therefore, the reverberation process on the accretion disk, warm absorbers, the 
BLR or the torus cannot explain the soft lags observed in Cen A. 

6 Summary 

Employing the DCF methods, we found evidence for soft lags on a timescale of days in the flux 
enhancement event from MJD 55075 to MJD 55150, between all energy bands’ X-ray variation in 2- 
4 keV, 4-10 keV, and 10-20 keV of Centaurus A in the X-ray light curve observed by the MAXI/GSC, 
and quantify the significance of these results by performing Monte-Carlo simulations known as the 
FR/RSS method. The confidence of existence of the day-scale soft lag is higher than 99%. Also, we 
find that the soft lag is mainly originated by flux variations on a timescale of ^ 20 days. 

If the observed X-rays are emitted by SSC in a single zone in the jet by a single injection 
with the parameters : B = 6.2 G, 7min = 300, and 5=1 given by Abdo et al. (2010), the day-scale 
soft lags is difficult to explain because of too short cooling time of highly energetic electrons. For 
interpreting the evaluated soft lags by the simple one zone SSC model, a multiple injection scenarios 
with a hard-to-soft evolution, or a weaker B field or/and a smaller Forentz factor of electrons are 
required. Alternatively, assuming that the variable X-rays originate in the disk photons which are 
up-scattered by electrons in a corona with a temperature of around 10 keV surrounding the inner part 
of the accretion disk, we show the possibility that the observed lags can be explained by the cooling 
of the Comptonizing corona. The plausibility of this scenario is supported by the consistency in the 
decay times and the break timescale of the power spectrum density. With the reverberation scenario, 
a hard lag is naturally expected, and therefore the observed soft lags cannot be explained. 
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Fig. 1. The five days binned long-term X-ray light curves of Centaurus A in the 2-4 keV (top panel), 4-10 keV (middle panel), and 10-20 keV (bottom panel) 
obtained by the maxi/GSC. 
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Fig. 2. CEDCFs and ZDCFs calculated between the 2-A keV and 4-10 keV, 4-10 keV and 10-20 keV, and 2-A keV and 10-20 keV for the long-term light 
curves. CFDCFs are shown as gray shadowed area, and ZDCFs are plotted by black circles with error bars. Dashed vertical line stands for ri^g = 0. For the 
definition of the horizontal error bars, see §3.3. 


18 










Probability 



Fig. 3. Cross correlation centroid distributions (CCCDs) ot CEDCF (solid line) and ZDCF (dashed iine) tor the long-term light curves obtained by the FR/RSS 
method. Derived soft iags and corresponding 1 a ranges are also shown above the histograms. Gray vertical lines are drawn at -nag = 0. 
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Fig. 4. Running CEDEFs for full length X-ray light curves between 2-4 keV and 4-10 keV, 4-10 keV and 10-20 keV, and 2-4 keV and 10-20 keV. The gray 
scale shows a correlation coefficient between different energy bands' light curves separated by 200days. The top panel is the 5 days binned light curve for 
reference. 
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Fig. 5. X-ray light curves of Centaurus A obtained by the MAXI/GSC. The flux data are plotted as the black points with error bars. The error bars present 1 o- 
errors derived from C-statistios (see Hiroi et al. 2011 for details). The bottom panel is the ratio of the 4-10 keV to the 2-A keV flux. 
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Fig. 6. CEDCFs and ZDCFs calculated between the 2-A keV and 4-10 keV, 4-10 keV and 10-20 keV, and 2-4 keV and 10-20 keV for the flux enhanced 
episode. CFDCFs are shown as gray shadowed area, and ZDCFs are plotted by black circles with error bars. Dashed vertical line stands for ri^g = 0. For 
the definition of the horizontal error bars, see §3.3. 
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Fig. 7. Cross correlation centroid distributions (CCCDs) of CEDCF (solid line) and ZDCF (dashed iine) for the flux enhancement event obtained by the FR/RSS 
method. Derived soft iags and corresponding f a ranges are also shown above the histograms. Gray vertical lines are drawn at -nag = 0. 
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